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ABSTRACT 

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has 
been observed. The complete ~ 3.3° x 3.3° imaging with SPIRE 250/350/500 yt/m and PACS 70/160 yt/m allows a deep investigation of embedded 
protostellar phases, probing of the dust emission from warm inner regions at 70 and 160yum to the bulk of the cold envelopes between 250 and 
500yL/m. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy 
distributions are derived to measure luminosities and envelope masses, and to place the protostars in an Menv - ^boi evolutionary diagram. The 
spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected 
protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with 
Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class YSOs newly revealed by Herschel. For the entire Aquila field, 
we find a total of ~ 45 - 60 Class YSOs discovered by Herschel. This confirms the global statistics of several hundred Class YSOs that should 
be found in the whole Gould Belt survey. 
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1. Introduction 

During the main accretion phase, protostars are deeply embed- 
ded in their collapsing envelopes and parent clouds. They are so 
embedded that they radiate mostly at long wavelengths, making 
their detection and study difl&cult from the ground (e.g. | Andre] 
|et al.|2000[|Di Francesco et al.|2007| ). Protostars, or young stellar 
objects (YSOs), in the solar neighborhood have been extensively 
surveyed, but a complete and unbiased census of all protostars 
in nearby molecular clouds is lacking. The census of embedded 
YSOs provided by IRAS and near-IR studies in the 1980s and 
1990s was far from complete even in the nearest clouds. Thanks 
to their high sensitivity and good spatial resolution in the mid- 
infrared, ISO and, more recently, Spitzer could perform more 
complete surveys in all major nearby star-forming regions (e.g. 



scale far-infrared images at an unprecedented spatial resolution 
and sensitivity. The Aquila Rift region has been chosen to be 
observed for the science demonstration phase of Herschel for 
this survey. 

Our 250/350/500 jim SPIRE and 70/160 jim PACS images 
of the Gould Belt provide the first access to the critical spec- 
tral range of the far-infrared to submillimeter regimes to cover 
the peak of the spectral energy distributions (SEDs) of the cold 
phase of star formation at a high enough spatial resolution to sep- 
arate individual objects. The Herschel surveys theref ore allow 
an unprecedented, unbiased census of starless cores ([Konyves 



|et al.||2010|), embedded protost ars (this work), and cloud struc- 
ture ( [Men'shchikov et al.|2010| ), down to the lowest column den 



% 



Nordh et al.l 19961 [Bontemps et al.|200lHKaas et al.|2QQ4l [Alien 



et al.||2QQ7| 



Evans et al.||2009| ). The population of the youngest 



protostars, the Class YSOs, can however not be properly sur- 
veyed solely in the near and mid-infrared. These youngest ob- 
jects remain weak or undetected shortward of ~ 20yum. 

The Herschel Gould Belt Survey ( [Andre et aLl|2010|) is a 
key prog ram of the ESA Herschel Space Observatory fP ilbratt 
[etal.[2010). It employs the SPIRE (GriffinetaL^OlO ) and PACS 
( [Poglitsch et al.[[2QTQ | instruments to do photometry in large- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 
** Figures [2}J3| are only available in electronic form at 



http ://w w w. aanda. org | 



sities ( Andre et al.[[2010| ). This survey yields the first accurate 
far-infrared photometry, hence good luminosity and mass esti- 
mates, for a comprehensive view of all early evolutionary stages. 



2. Observations 

The observations were performed in the parallel mode of 
Herschel with a scanning speed of 60'7sec, which allows photo- 
metric imaging with SPIRE at 250, 350, and 500 yum and PACS 
at 70 and 160yum. Two cross-linked scan maps were performed 
for a final coverage of ~ 3.3° x 3.3° (see Fig.[T]). 

The SPIRE data were reduced using HIPE version 2.0 and 
modified pipeline scripts; see GrifiS n et al.[ ( [2Q1Q|) for the in- 
orbit performance and scientific capabilities, and |Swinyard et al. 
( |2010| ) for calibration methods and accuracy. A median baseline 
was applied to the maps for each scan leg, and the naive map- 
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Fig. 1. Visual extinction map towards the whole Aquila Rift / 
Serpens region derived by us from 2MASS data (see details 
in Schneider et al._2Q10K The spatial resolution is T FWHM. 
The dashed blue rectangle indicates the Herschel coverage. It 
comprises the bright Hll region W40 (white circle), the Serpens 
South cluster (white star), and the Hll region Sh2-62, associated 
with the young star MWC297. The Aquila Rift corresponds to 
the large elongated structure from the northeast to the southwest. 
Emisson from the Galactic plane is seen in the southeast comer. 



per was used as a mapmaking algorithm. The RACS data were 
reduced in HIRE 3.0. We used an updated version of the cali- 
bration files following the most recent prescriptions of the RACS 
ICC (see |Konyves et al.|2QlQ| for details). Multiresolution me- 
dian and second-order deglitching, as well as a high-pass filter- 
ing over the full scan leg length, were applied. The final RACS 
maps were created using the photRroject task, which performs 
simple projection of the data cube on the map grid. 

The resulting RACS maps are displayed in Fig.|2] Owing to 
the rapid mapping speed, the resulting point spread functions 
(RSFs) are elongated in the scan directions leading to cross- 
like shapes of the RSFs with expected sizes of 5. 9'' x 12.2'' at 
70yum and 11.6''xl5.7'' at 160 yum. The resulting rms in these 
maps ranges from 50 to 1000 mJy/beam at TOyum and from 120 
to 2200 mJy/beam at 160yum, depending on the level of back- 
ground in the map. It is in the W40/Sh2-64 Hll region that the 
background level is the highest. 



On the other hand, the most active and main extinction fea- 
ture in the 2MASS extinction map is associated with the HII 
region W40/Sh2-64, which has so far been considered to be 
at a distance ranging from 100 and 7 00 pc depending on au- 
thor dSmith et aLl[l985l |Vallee|[T987 ^ and references therein). 
These distance estimates are mostly based on kinematical dis- 
tances that have large uncertainties. W40 could therefore beat 
the same distance as Serpens. Recently, Gutermuth et al. (20081 
reported Spitzer observations of an embedded cluster, referred 
to as Serpens South, in the Aquila Rift region. This cluster is 
located very close in projection on the sky to W40 (see Fig. \\\ 
and thus seems to be part of the W40 region. Guterm uth et al. 
(2008 ) proposed that the Serpens South cluster should be part 
of Serpens since it has the same velocity (6 km/s). The molecu- 
lar cloud associat ed with W40 and traced by CII recombination 
lines and CO (Zeil ik & Lada|1978] ) has a velocity ranging from 
4.5 to 6.5 km/s, which is also roughly the same as Serpens. More 
recent N2H^ observations of the entire W40/Serpens South re- 
gion confirm similar velocities in the whole region with velocity 
diff'erences of only ~ 2 km/s (Maury et al. in prep). It is therefore 
more straightforward to consider that the W40 region is a sin- 
gle complex at the same distance as Serpens. This distance also 
suits the MWC297 / Sh2-62 region since the young IOMq star 



3. Overview and distance of the Aquila Rift complex 

The Aquila Rift is a coherent, 5° long feature above the Galactic 
plane at 1=28°, clearly visible on an extinction map derived 
from the reddening of stars in 2MASS (Fig. [T]). A distance 
of 225 ± 55 pc has been derived for this extinction wall, us- 
ing spectro-p hotometric studies of the optically visible stars 
( Straizys et al. 2003). This distance is very similar to the usu- 
ally adopted distance of 260 ± 37 pc for the S erpens star-forming 
regioiQ located only 3° north ( Straizys et al._1996) . 



^ Note that a larger distance of 415 + 25 pc has been recently claimed 
for Serpens Main based on a VLBA parallax of EC95, a young AeBe 
star embedded in Serpens Main (|Dzib et al. 2010|). 



MWC297 itself has an accepted distance of 250 pc ( [Drew et aL] 
|1997|). It is finally worth noting that the visual extinction map 
by [Cambresy (1999) derived from optical star counts and only 
tracing the first layer of the extinction wall has exactly the same 
global aspect as the 2MASS extinction map of Fig. [T] suggesting 
that both Serpens Main and the W40 / Aquila Rift / MWC297 re- 
gion are associated with this extinction wall at 260 pc. We thus 
adopt the distance of 260 pc for the entire region in the follow- 
ing. 

The 2MASS extinction map and the Herschel images (see 
Appendix for RACS images and Ko nyves et al.|2010| for SRIRE 
images) clearly show a massive cloud associated with W40. This 
cloud corresponds to G28.74-h3.52 in ( [Zeilik & Lada|1978| ) and 
has a mass of 1.1 x 10^ M© (derived from our 2MASS extinction 
map). The cloud associated with MWC297 is less massive (4.1 x 
10^ Mq), and we obtain a total mass of 3. 1x10"^ M© for the whole 
area covered by Herschel. 



4. Results and analysis 

4.1. Source detection and identification of protostars 

A systematic source det ection was performed on all 5 Herschel 
bands using getsources ( [Men'shchikov et aT][2010| ). This code 
uses a method based on a multiscale decomposition of the im- 
ages to disentangle the emission of a population of spatially co- 
herent sources in an optimized way in all bands simultaneously. 
We built a sample of the best candidate protostars for the whole 
field. These sources are clearly detected in all Herschel bands 
(high significance level), and we require a detection at the short- 
est Herschel wavelength, 70yum (or 24yum when Spitzer data 
were available), to distinguish YSOs from starless cores. Since 
the 24 and 70yum emission should only trace warm dust from the 
inner regions of the YSO envelopes, these sources can be safely 
interpreted as protostars. The 70yum fluxes have even been re- 
cently recognized as a very good tracer of protostellar luminosi- 
ties ( Dunham et al.|20Q8 ). On the other hand, in the RDR region 
of W40 some extended emission from warm dust at the Hll re- 
gion interface could contaminate this YSO detection criterium. 
To avoid too stringent a contamination from this extended emis- 
sion, we selected only sources with an FWHM size smaller than 
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Fig. 4. Spectral energy distributions of a newly discovered bright Class (left panel) and of a weaker Class object (right panel), 
which was previously detected with Spitzer in|Gutermuth et al.|(|2008|). 



4C at 70yum. Also, we had to make a source detection using a 
large pixel size of 6'\ which is good enough for starless cores 
mostly detected in the SPIRE bands but not perfect to sample 
the spatial resolution at 70yum and properly disentangle possi- 
ble multiple protostellar sources. A more precise detection could 
only be achieved in a reduced area in the Serpens South region 
(see Sect. [44]). 

A large number of compact sources are clearly seen in the 
70yum map down to the sensitivity limit of the survey. In the 
whole Aquila field, 201 YSOs were detected with getsources. 
The best achieved rms (50 mJy/beam) in the Aquila TOyum map 
in the lowest background regions corresponds to a 5 cr detec- 
tion level in terms of protostar luminosity of 0.05 L© using the 
[Dun ham et al.^ (2008) relationship. In contrast, in the highest 
background regions, the 5 cr detection level is then as high as 
1.0 Lq. To account for the variable background level in Aquila, 
we performed simulations to evaluate the final completeness 
level of the YSO det ection and obtained a 90 % completeness 
level of -0.2 L© (see Konyves et al."2010), which is compatible 
with the above rough estimates using Dunham et al. (2008). 




3 30 40 20- 
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4.2. Spatial distribution of the protostars 

We plotted in Fig.[3]the spatial distribution of the Herschel sam- 
ple of 201 YSOs overlaid on the map of the dust temperature 
derived from a simple graybody fit of the Herschel data (see de- 
tails in Konyves et al.pOlO] ). It is clear that the W40 region cor- 
responds to the most active star-forming region in the Herschel 
coverage with 90 % of the detected protostars. A second, much 
less rich, site corresponds to MWC297 with 8 % of the proto- 
stars, and another site to the east of W40 can be tentatively iden- 
tified with very few candidate protostars. 

4.3. Basic properties of the protostars: Menv and Lboi 

For each source an SED was built using the 5 bands of Herschel, 
as well as Spitzer photometry (Gutermu th et al.| [2008) and 
MAMBO 1.2mm data (Maury et al. in prep) when available. 
These SEDs were systematically fitted using graybody functions 
to derive Menv in a systematic way, while the basic properties 
Lboi, ^^ubmm' ^^^ -^boi wcrc obtained by simple integrations of 
the SEDs. Two representative SEDs are displayed in Fig. [4] with 
a newly discovered Class object and a weaker Class s ource, 
which has a Spitzer counterpart in [Gutermuth et al. ( 2008| ). 



Fig. 5. Herschel SPIRE 350yum image of the Serpens South re- 
gion with the distribution of Herschel candidate protostars (blue 
circles) from the whole field extraction, of the 7 newly discov- 
ered Class protostars (green circle s), and of the Spitzer YSOs 
(red crosses; [Gutermuth et al.|2008] ). 



4.4. A close-up view of the Serpens South region 

To go one step further in the identification and characterization 
of the Herschel protostars, we performed a more detailed analy- 
sis of the sources in a small area around the Serpens South clus- 
ter. In this area, we made a dedicated getsources source extrac- 
tion using a smaller size pixel of y\ and we could compare t hese 
first results with the Spitzer protostar population by Guter muth[ 
[et al., ( ,2008) . We used getsources on 8 bands from 8 to 1200//m 
by adding the 8 and 24yum Spitzer and the 1.2mm MAMBO data 
to the 5 Herschel bands. 

A synthesized view of these first results based on this novel 
panchromatic analysis of infrared to millimeter range data for 
this area is given in Fig. [5] It shows the distribution of Herschel 
protostars compared to the Spitzer sources. The first analysis of 
this field indicates that even in a highly clustered region like 
Serpens South, a significant population of protostars were found 
to be missing by pure near and mid-infrared imaging with as 
many as 7 newly detected Class objects in this field. We also 
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note that the Spitzer protostars (most of these not detected with 
Herschel) probably correspond to evolved or low-luminosity 
Class I objects. 



5. Global view of the protostellar population in 
Aquila 



Using the basic properties derived in Sect. |4.3[ we can draw the 
first picture of the property space Herschel is going to cover 
thanks to its unprecedentedly sensitive and high spatial resolu- 
tion in the far-infrared. 

In Fig. [6] we plotted the location of the 201 Herschel YSOs 
obtained in the entire field in a Mgnv - ^boi evolutionary diagram 
used to compare observed properties with theoretical evolution- 
ary models or tracks. The displayed tracks represent the expected 
evolution of protostars of masses 0.2, 0.6, 2.0, and 8.0 M© from 
the earliest times of accretion (upper left part of the diagram) 
to the time of 50 % mass accreted (conceptual limit between 
Clas s and Class I YSOs), and the time for 90 % mass accreted 
(see Bontemps et al. 1996[ Saraceno et al. 1996[ Andre et al. 
2000[ [Andre et al."2008'). In this plot, we distinguished objects 



with an ^^^'^V^bof'''''' ^^^^^^ ^^^^ ^'^^ which could be safely 
recognized as Class objects, from YSOs with J-^^^^J JJ^^^^^^ 
lower than 0.01 which are proposed to be Class I sources. The in- 
termediate objects with 0.03 > ^^ubmrn/^bof^^^ > ^-^^ should be 
seen as objects with an uncertain classification. A forthcoming 
analysis will resolve their nature by building complete SEDs in- 
cluding Spitzer data for a large part of the Aquila field. So far we 
could safely classify objects only in the reduced area of Serpens 
South (Sect. |4.4|). In this subfield, we verified that objects with 



submm^ ^bol 



7A 



70-5 



0.03 and T^-^^^ 



< 27 K using the reduced 



(only the 5 Herschel bands from 70 to 500 yum) SED coverage 
are indeed all found to be Class objects based on the full cov- 
erage from Syum to 1.2 mm. We see that the obtained location 
of Class and Class I YSOs is compatible with the 50 % mass 
accreted limit. Imposing Class objects to have to be above this 
limit (dashed line in Fig. [6]), we finally found between 45 (for 
J70-500 ^ 27 K) and 60 (^^bmm/^boi"^^^ > 0-^^) ^^^^^ ^ ^^J^^^s 
in the entire field of Aquila. 

In conclusion, even if the precise locations of the Herschel 
protostars in this diagram are seen as a preliminary result and 
will be updated with a more complete analysis and source detec- 
tion, our early results clearly indicate that Herschel is a power- 
ful tool for probing the virtually unexplored area of the physical 
properties of the earliest stages of protostellar evolution. 
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Fig. 2. PACS 70 /im (left) an d 160)um (right) images of the Aquila field. See details about data redu ction and map making in Sect.|2| 
and in Konyves et al. ( 2010| ). The corresponding SPIRE 250, 350, and 500 yum images are shown in Konyves et al. ( 2010| ). 
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Fig. 3. Distribution map of the 201 Herschel YSOs selected in Sect. |4.1| over-plotted on the map of dust temperature. The dust 
temperature map was derived from graybody fits to the Herschel data (see details in |Konyves et al.|2010| ). 



